arXiv: 1507.03200vl [quant-ph] 12Jul2015 


Submitted manuscript No. 

(will be inserted by the editor) 


Duality quantum computer and the efficient 
quantum simulations 


Shi-Jie Wei • Gui-Lu Long 


July 12, 2015 


Abstract In this paper, we firstly briefly review the duality quantum com¬ 
puter. Distinctly, the generalized quantum gates, the basic evolution oper¬ 
ators in a duality quantum computer are no longer unitary, and they can 
be expressed in terms of linear combinations of unitary operators. All lin¬ 
ear bounded operators can be realized in a duality quantum computer, and 
unitary operators are just the extreme points of the set of generalized quan¬ 
tum gates. A d-slits duality quantum computer can be realized in an ordinary 
quantum computer with an additional qudit using the duality quantum com¬ 
puting mode. Duality quantum computer provides flexibility and clear physical 
picture in designing quantum algorithms, serving as a useful bridge between 
quantum and classical algorithms. In this review, we will show that duality 
quantum computer can simulate quantum systems more efficiently than or¬ 
dinary quantum computers by providing descriptions of the recent efficient 
quantum simulation algorithms of Childs et al [Quantum Information & Com¬ 
putation, 12(11-12): 901-924 (2012)] for the fast simulation of quantum sys¬ 
tems with a sparse Hamiltonian, and the quantum simulation algorithm by 
Berry et al [Phys. Rev. Lett. 114 , 090502 (2015)], which provides exponential 
improvement in precision for simulating systems with a sparse Hamiltonian. 
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1 Introduction 

One of us. Long, came to know Dr. Brandt first through his important works 
in quantum information [Tl[2jl3l|4l[5l[6j[7l|8l[9] and later his role as editor-in-chief 
of the journal Quantum Information Processing(QIP). Long proposed a new 
type of quantum computer in 2002 nni, which employed the wave-particle du¬ 
ality principle to quantum information processing. His acquaintance with QIP 
began in 2006 through the work of Stan Gudder who established the math¬ 
ematical theory of duality quantum computer which was accompanied 
by a different mathematical description of duality quantum computer in the 
density matrix formalism m- The development of duality quantum computer 
owes a great deal to QIP first in the term of Dr. David Cory as editor-in-chief, 
and then the term of Dr. H. E. Brandt as the editor-in-chief. For example, the 
zero-wave function paradox was pointed out firstly by Gudder m, and two 
possible solutions were given in Refs. m and M- Long has actively partici¬ 
pated in the work of QIP as a reviewer when Dr. Brandt was the editor-in-chief, 
and as a member of the editorial board from 2014. At this special occasion, 
it is our great honor to present a survey of the duality quantum computer in 
this special issue dedicated to the memory of Dr. Howard E. Brandt. 

As is well-known, a moving quantum object passing through a double¬ 
slit behaves like both waves and particles. The duality computer, or duality 
quantum computer exploits the wave-particle duality of quantum systems m- 
It has been proven that a moving n-qubit duality computer passing through 
a d-slits can be perfectly simulated by an ordinary quantum computer with 
n-qubit and an additional d levels qudit [T51IT61IT7] . So we do not need to build 
a moving quantum computer device which is very difficult to realize. This 
also indicates that we can perform duality quantum computing in an ordinary 
quantum computer, in the so-called duality quantum computing mode [Bin]- 
There have been intensive interests in the theory of duality computer in recent 
years . and 

experimental studies have also been reported mm- 

This article is organized as follows. In section[21 we briefly describe the gen¬ 
eralized quantum gate, the divider and combiner operations. Section [3] reviews 
the duality quantum computing mode, which enables the implementation of 
duality quantum computing in an ordinary quantum computer. In section|4l we 
outline the main results of mathematical theory of duality quantum computer. 
In section [21 we give the duality quantum computer description of the work of 
Ghilds et al [33] which simulates a quantum system with sparse Hamiltonian 
efficiently. In section [51 we give the duality quantum computer description of 
the work of Berry et al m which simulates a quantum system having a sparse 






Duality Quantum Computing 


3 


0 


1 


2 




Detector 


Fig. 1 An illustrated picture for a three-slits duality quantum computer. The input is from 
the 0-th slit, and the output of duality quantum computing is taken from only Oth-slit on 
the right wall m- 


Hamiltonian with exponential improvement in the precision. In section [71 we 
give a brief summary. 


2 Duality quantum computer, Divider, Combiner operations 

A duality quantum computer is a moving quantum computer passing through 
a d-slitsflU]. In Fig. [I] we give an illustration for a duality quantum computer 
with 3-slits m- The quantum wave starts from the 0-th slit in leftmost wall, 
and then goes to the middle screen with three slits (this is the divider op¬ 
eration). Between the middle screen and the rightmost screen, some unitary 
operations are performed on the sub-waves at different slits. They are then 
collected at the 0-slit in the rightmost screen, and this is the quantum wave 
combiner operation. The result is then read-out by the detector placed at the 
0-slit on the right wall. 

In a duality quantum computer, the two new operations are the quantum 
wave divider (QWD) operation and quantum wave combiner (QWC) operation 
[10) . The QWD operation divides the wave function into many identical parts 
with reduced amplitudes. The physical picture of this division is simple and 
natural : a quantum system passing through a d-slits with its wave function 
being divided into d sub-waves. Each of the sub-waves has the same internal 
wave function and are different only in the center of mass motion characterized 
by the slit. Conversely, the combiner operation adds up all the sub-waves into 
one wave function. It should be noted that one divides the wave function of 
the same quantum system into many parts in quantum divider, whereas in 
quantum cloning one copies the state of one quantum system onto another 
quantum system P^I3h) (which may also holds true for classical systems [32) ■ 
So, the division operation does not violate the non-cloning theorem. 
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Considering a quantum wave divider corresponding to a quantum system 
passing through a d-slits. Writing the direct sum of Hilbert space as 

the form of H® where Hi = H, i G {0, ••• ,d — 1}. The divider structure 
characteristics which describes the properties of a quantum wave divider can 
be denoted as {pi, * = 0, • • • , d — 1} where each pi is a complex number with a 
module less than 1 and satisfy X]i=i biP = 1- The divider operator Dm which 
maps H ^ H® is defined by 

Dm-ip = <BilliPii^i). ( 1 ) 

where ipi = ip^ i G {0^ ■ ■ ■ ,d—1}. This is the most general form of the divider 
operator, and it describes a general multi-slits. In a special case, the multi-slits 
are d identical slits, then pi = \/TJd. 

The corresponding combiner operation Cm can be defined as follows 


d-l 

0 • * * © ^d-i) = ^ ( 2 ) 

where q = {gor ’ ’ idd-i} denotes the combiner structure that describes the 
properties of a quantum wave combiner. Each qi is a complex number that 
satisfy the module less than 1, and = 1. In the case of qi = \/TJd, the 

combiner structure is uniform. 

Now, we consider the uniform divider and combiner structures which cor¬ 
respond to Pi = \/T/d and qi = \/TJd, respectively. In this case, the combined 
operations of divider and combiner leave the state unchanged. The process can 
be described as follows 

CmDmIp = Cm ©fcl tpi = '^ 1 /d tjj = -Ip. (3) 

i 

If the divider structure and combiner structure satisfy certain relation, this 
property also holds. The details will be given in the next section. 

It will be shown later in the next section that the divider and combiner 
structure Dm and Cm can be expressed by a column or a row of elements of 
a unitary matrix respectively. For duality quantum gates with the form of Lr 
in Eq. (|6]) , the relation of the two unitary matrices makes the structures of Cm 
and Dm adjoint of each other. 


3 Duality Quantum Computing Mode in a Quantum Computer 

The most general form of duality quantum gates has been given in Refs. [161 
in]. For the convenience of readers, we use the expressions from duality quan¬ 
tum computing mode nsumna. Compared to ordinary quantum computer 
where only unitary operators are allowed, the duality quantum computer offers 
an additional capability in information processing: one can perform different 
gate operations on the sub-wave functions at different slits m- This is called 
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the duality parallelism, and it enables the duality quantum computer to per¬ 
form non-unitary gate operations. The generalized quantum gate, or duality 
gate is defined as follows 

d-l 

L,=J2c^U,, (4) 

where Ui is unitary and q is a complex number and satisfies 

d-l 

( 5 ) 

i=0 

The duality quantum gate is called real duality gate or real generalized quan¬ 
tum gate when it is restricted to positive real c^. In this case, Ci is denoted by 
Ti, and they are constrained by the condition of = 1- The real duality 

gate is denoted as Lr- So, the real duality quantum gate can be rewritten as 

d-l 

Lr = Y,nU,. ( 6 ) 

i=0 

This corresponds to a physical picture of an asymmetric d-slits, and is the 
probability that the duality computer system passes through the i-ih slit. 

According to the definition of duality quantum gates, they are generally 
non-unitary. It naturally provides the capability to perform non-unitary evo¬ 
lutions. For instance, dynamic evolutions in open quantum systems should be 
simulated in such machines. More interestingly, it is an important issue to 
study the computing capabilities of duality quantum computing. An impor¬ 
tant step toward this direction is that Wang, Du and Dou [19] proposed an 
theorem which limits what can not be a duality gate in a Hilbert space with 
infinite degrees of freedom. 

The divider operation can be expressed by a general unitary operation V 
and the combiner operation can be expressed by another general unitary op¬ 
eration W. The two unitary operations are implemented on an auxiliary qudit 
which represents a d-slits. The quantum circuit of duality quantum computer 
is shown in Fig. [5| There are d controlled operations between the operations 
V and W. The d energy levels of the qudit represent the d-slits. 

We divide the duality computing processing into four steps to reveal the 
computing theory in a quantum computer. 

Step one: Preparing the initial quantum system |'F)|0) where |<F) is the 
initial state. Then performing the divider operation by implement the V on 
the auxiliary qudit, and this operation transforms the initial state into 

d-l 

\'l')\0)^\dr)V\0)=Y,V.omi)- (7) 

i=0 

Vio is the first column element of the unitary matrix V representing the co¬ 
efficient in each slit. V^o = Pi-, denotes the divider structure. Note that Vio 
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Fig. 2 The quantum circuit of duality computing in a quantum computer. denotes the 
initial state of duality quantum computer with |0) as the controlling auxiliary qudit. The 
squares represent unitary operations and the circles represent the state of the controlling 
qudit, Unitary operations Uq, Ui ■ • ■, Ud-i are activated only when the qudit holds the 
respective values indicated in circles m- 


is a complex number with |Vio| < 1 and satisfies = 1- So V is a 

generalized quantum division operation. 

Step two: We perform the qudit controlled operations Uq, Ui - ■ ■, Ud-i on 
the target state jlf') which leads to the following transformation 


d-l 

( 8 ) 

i=0 

This corresponds to the physical picture that implements unitary operations 
simultaneously on the sub-waves at different slits. 

Step three: We combine the wave functions by performing the unitary 
operation W. The following state is obtained, 

Y,V^oU,mw\i). (9) 


Step four: Detecting the final wave function when the qudit is in state |0) 
by placing a detector at slit 0 as shown in Fig. [2] The wave function becomes 

Y.Wo^VMW)\0) = ^(VFo.V.o)C/^|‘F)|0). (10) 

i i 

It should be noted that TV is a generalized quantum combiner operation and 
qi = Woi, which is the combiner structure in Eq. Hence Ci = W^iVio is 
the coefficient in the generalized duality gate defined in Eq. ©• Now, we have 
successfully realized the duality quantum computing in an ordinary quantum 
computer. 

Considering a special case that W = V\ the coefficients Cj = = 

|VioP satisfy 




where is defined in Eq. m corresponding to the real duality gate Lr- 
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Generally speaking, Ci = WoiVio is a complex number. The sum of q’s can 
be denoted as 

= (W)oo. (12) 

i i 

The value of (ITy)oo is just an element of a unitary matrix WV, and naturally 
has the constraint |(lTy)oo| < 1- Hence the most general form of duality gates 
allowable by the principles of quantum mechanics is the form of ([S]) 

In a recent study, Zhang et al [35] has proved that it is realizable and 
necessary to decompose a generalized quantum gate in terms of two unitary 
operators V and W in Eq. (IT2]l if and only if the coefficients satisfy l^il < 1- 
Obtaining the explicit form of the decomposition is a crucial step in duality 
quantum algorithm design and related studies. 


4 Mathematical theory of duality quantum computer 

The mathematical theory of duality quantum computer has been the subject 
of many recent studies [TUHTTlfT^lfTBlfTTlITSlIT^ll 2 0 II 2T1I33] . Here we briefly review 
the mathematical description in duality quantum computing. In this case, the 
mathematical theory of the divider and combiner operations are restricted to 
a real structure, namely each pi{i = 0,1, • • • , d — 1) is real and positive, and 
the uniform structure is also a special case of the combiner structure. The 
following results are from Ref. m and we label the corresponding lemma, 
theorem and corollaries by a letter G, and the corresponding operators are 
labeled with a subscript p. 

Here are the properties of generalized quantum gates and related operators 

[Si- 

Defining the set of generalized quantum gates on H as G{H) which turns 
out to be a convex set. Then we have 

Theorem G 4. 1 The identity Ih is an extreme point of G{H), where Ih 
is the identity operator on H. 

Any unitary operator is in Gi'H) and In € G{T~L). The identity Ih is an 
extreme point of G{H) which indicates 'Y^^PiUi = Ih if and only if Ui = Ih 
for all i. 

Corollary G 4. 2 The extreme points of Gi'H) are precisely the unitary 
operators in H. 

We can conclude from Theorem G 4. 1 and Corollary G 4. 2 that 
the ordinary quantum computer is included in the duality quantum computer. 

Denoting B{'H) by the set of bounded linear operators on H and let 
be the positive cone generated by GiH). That is 

M+G{H) = {aA-. A&GiH),a>Q}. (13) 

Theorem G 4. 3 If dim'H < oo, then i3('H) = R+t/)?^). 

This theorem shows us that the duality quantum computer is able to sim¬ 
ulate any operator in a Hilbert space H if dim'H < oo. 
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It should be pointed out that these lemmas, corollary and theorems also 
hold for divider and combiner with a general complex structure m- One 
limitation has been given explicitly by Wang, Du and Dou [19] that what 
can not be a generalized quantum gate when the dimension is infinite. It is 
an interesting direction to study the computing ability of duality quantum 
computer in terms of this theorem. 


5 Description of Childs-Wiebe Algorithm for Simulating 
Hamiltonians in a Duality Quantum Computer 

Simulating physics with quantum computers is the original motivation of 
Richard Feynman to propose the idea of quantum computer m- Benioff has 
constructed a microscopic quantum mechanical model of computers as repre¬ 
sented by Turing machines [39]. Quantum simulation is apparently unrealis¬ 
tic using classical computers, but quantum computers are naturally suited to 
this task. Simulating the time evolution of quantum systems or the dynamics 
of quantum systems is a major potential application of quantum computers. 
Quantum computers accelerate the integer factorization problem exponentially 
through the use of Shor algorithm [40] , and the unsorted database search prob¬ 
lem in a square-root manner through the Grover’s algorithm [41] (see also the 
improved quantum search algorithms with certainty [421143] 1. Quantum com¬ 
puters can simulate quantum systems exponentially fast [3B]. Lloyd proposed 
the original approach to quantum simulation of time-independent local Hamil¬ 
tonians based on product formulas m which attracted many attentions [13 
I46ll47ll4^ ?]. However, in this formalism, high-order approximations lead to 
sharply increased algorithmic complexity, the performance of simulation algo¬ 
rithms based on product formulas is limited [33] . For instance, the Lie-Trotter- 
Suzuki formulas, which is high-order product formulas, yields a new efficient 
approach to approximate the time evolution using a product of unitary oper¬ 
ations whose length scales exponentially with the order of the formula [51] . In 
contrast, classical methods based on multi-product formulas require a sum of 
unitary operations only in polynomially scales to achieve the same accuracy 
[5^ . However, due to the unclosed property of unitary operations under ad¬ 
dition, these classical methods cannot be directly implemented on a quantum 
computer. 

The duality quantum computer can be used as a bridge to transform clas¬ 
sical algorithms in to quantum computing algorithms. Duality parallelism in 
the duality quantum computer enables us to perform the non-unitary oper¬ 
ations. Moreover,duality quantum gate has the form Lc = YltZo ^iUi- This 
is the linear combinations of unitary operations. Duality quantum computer 
is naturally suitable for the simulation algorithms of Hamiltonians based on 
multi-product formulas. 

Childs and Wiebe proposed a new approach to simulate Hamiltonian dy¬ 
namics based on implementing linear combinations of unitary operations [331 
[55] . The resulting algorithm has superior performance to existing simulation 
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algorithms based on product formulas and is optimal among a large class of 
methods. Their main results are as follows 

Theorem 1 Let the system Hamiltonian be H = where each Hj C 


is Hermitian and satisfies || Hj ||< h for a given constant h. Then the 
Hamiltonian evolution can he simulated on a quantum computer with 


failure probability and error at most e as a product of linear combinations of 
unitary operators. In the limit of large to, ht, 1/e, this simulation uses 


6 ( TO2htei ®Vi°g("^'»*/E) 


(14) 


elementary operations and exponentials of the HjS. 

Considering this simulation algorithm is based on implementing linear com¬ 
binations of unitary operations, it can be implemented by duality quantum 
computer. Now, we give the duality quantum computer description of this 
simulation algorithm. 

The evolution operator U (t) satisfies the Schrodinger equation 


ijUit) = HU{t) 


(15) 


and time evolution operator can be formally expressed as U{t) = 

The Lie-Trotter-Suzuki formulas approximate time evolution operator U 
for H = X)yLi Hj as a product of the form 


m 



These formulas can be defined for any integer x > 0 by [331I5T] 


1 


m 


5i(t) = 



'S'xW = ('S'x-i(sx-it))^ - 4 sx-i ]0 > ( 16 ) 


where Sy-i = (4 — 4^/*^^^ ^)) ^ for any integer x > 1. This choice of S;^_i is 


made to ensure that the Taylor series of matches that of e to 0(t^^+^). 


With the values of x large enough and the values of t small enough, the 
approximation of U(t) can reach arbitrary accuracy. 

Childs et al have simulated U{t) using r iterations of Mk^k{t/r) for some 
sufficiently large r m- 




( 17 ) 
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where £q {q € {1,2,..., fc+1}) represent distinct natural numbers and G 

1 (Cl,..., Ck+i G i?). In [33] , the £q and Ci are defined as 


£q 


g, if g < k, 

g7(fc+i)^ if g = /c + 1, 


and 





2T{fc + l) 

g2',(fc + l) 

1 ) 


2 ) 


if g < fc, 
if g = fc + 1. 


The formula is accurate to order, namely, 


||Mfc,fc(A)-[/(A)|| G 0(^4'=+!). 


(18) 


(19) 


( 20 ) 


The basic idea of this simulation algorithm is that dividing evolution time 
t into r segments and approximating each time evolution operator segment 
U{t/r) = by a sum of multi-product formula, namely, 

fc-i-i 

U{t/r) « MkAt/r) =Y,CqS^{t/£qrY-. (21) 

9=1 

Now, we give a duality quantum computer description of the implemen¬ 
tation of this simulation algorithm of time evolution. The quantum circuit is 
the same as Fig. |31 According to (ITHl) . Skit/iqrY” is an unitary operation. Let 
Ui = Sk{t/iqrY'^ and Ci = Cq, Mk,k{t/r) can be rewritten as 

k 

Mk,k{tlr) = ^ CiUi ( 22 ) 

2 = 0 


where Ui is an unitary operation. 

It is obvious that Mk^kit/r) is a duality quantum gate. The QWD is sim¬ 
ulated by the unitary operation V and the QWC is simulates by unitary op¬ 
eration VF on a qudit. The auxiliary qudit controlled operations is Ui. The 
matrix element Fo of the unitary matrix V and the matrix element Wio of the 
unitary matrix W satisfy: 


F.o = 


VFo,. = 


Vi,o 


^/E 


Vi,o 
Wo,i 


2 ’ 


As defined in (HU, Ci is the product of two unitary matrix elements: 


(23) 

(24) 


Ci — — 


^ 2 , 0 ^ 0,2 

\/E* 


(25) 
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The sum of q’s is 

= (W)oo. (26) 

i i 

In the special case W = V\ the simulation algorithm has the maximum 
success probability. The expression of Vi^o and 1Toy can be simplified into the 
form: 


14.0 = Wo,^ 


TETR' 


(27) 


After implementing the QWD operation, the auxiliary qudit controlled 
operations and the QWC operation, detecting the final wave function when 
the auxiliary qudit is in state |0). The initial state |!f')|0) has been transformed 
into 


k 

Y.Wo^V,oU^\<^')\0)=J2{Wo^V,o)U^\<^')\0)=Y,c^U,\^)\0). (28) 

i i 2=0 


The approximated evolution operator is implemented success¬ 

fully by the duality quantum computer. Implementing r segments of Mk^k{t/r), 
we can get the approximation of U{t) by Mk,k{t/rY■ Thus, this algorithm is 
clearly realized by the duality quantum computer in straightforward way. The 
essential idea of this algorithm is an iterated approximation, with each con¬ 
trolled Ui adding an additional high order approximation to the evolution 
operator. 


6 Description of Berry-Childs quantum algorithm with exponential 
improved precision for a sparse Hamiltonian system 

Berry and Childs provided a quantum algorithm for simulating Hamiltonian 
dynamics by approximating the truncated Taylor series of the evolution opera¬ 
tor on a quantum computer |34| . This method is based on linear combinations 
of unitary operations and it can simulate the time evolution of a class of phys¬ 
ical systems. The performance of this algorithm has exponential improvement 
over previous approaches in precision. 

Hamiltonian can be decomposed into a linear combinations of unitary op¬ 
erations: 

L 

H = (29) 

i=i 

Dividing the evolution time into r segments of length t/r. The time evolution 
operator of each segment can be approximated as 

Ur ■■= exp(-fiJt/r) Ri ^ ^ 
fe =0 
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where the Taylor series is accurate to order K. Substituting the Hamiltonian 
in terms of a sum of Hi into dsni), we can rewrite the truncated Taylor series 
as [34] 


Ur^U = f2 E ^ (31) 

fc=0 = l 

For convenience, we can set each > 0. Considering Hi is an unitary 
operation, we can conclude that the approximation H is a linear combinations 
of unitary operations. The expression has a quantum duality gate form. The 
truncated Taylor series index can be defined as [34] 

J := {(fc,£i,...,4) e J : fc < if}. (32) 

Then, the expression of U can be simplified as 

4 = (33) 

jeJ 

where ■= [i^ / rf / k']ai^ ■ ■ ■ ai^ ■= {-if ■ ■ ■ Hi^. 

It should be noted that U is not normalized. 

We define the normalization constant as s = According to ([S]), 

Lr = U/s is a. quantum duality gate. We let = j3i/s = Ui^ then it comes 

back to the duality quantum gate form in ([S]), 


d-l 

Lr = E 

i=0 


(34) 


where i G {0,1, 2,..., d — 1}. 

To give the duality quantum computer description, we need to realize the 
following processing 


|'F)|0)^|i')C/|0). (35) 

In Fig. [21 we give an illustration for our method to perform the algorithm 
in the form of quantum circuit. The unitary operation C/q corresponds to the 
decomposing form of Hamiltonians: H = OitHi and the quantum circuit 

In Fig. [3] implements Ur = ex.p{—iHt/r) ps J2k=o ■ 

The implementation of operation U need an L level auxiliary qudit |0 )l 
and K auxiliary qubits |0)^ which correspond to implementation of two QWD 
operations and QWC operations. Actually, the equation of (1311) indicates that 
we need summarize twice to realize the right side of this equation. We express 
the initial state as |>F)|0)^|0 )l 
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Fig. 3 Duality computing in a quantum computer. \^) is the initial state of duality quantum 
computer and there are /Cnumbers of |0) auxiliary controlling qubit and an L level auxiliary 
qudit|0)£, auxiliary controlling qudit. Unitary operations Uq are activated only when the 
qubit and qudit holds the respective values indicated in circles. The each unitary operation 
Uq is composed of Hi, H 2 , •.., Hl_i, 


Firstly, we transform the |0)^ part of the initial state into the normalized 
state using the QWD operation. We have 

^ I - 

(36) 

k=0 


We let / = J2o |t and define the QWD as which can be expressed as a 
2^ X 2^ matrix. The elements of the matrix satisfy 


^i,o — 



(37) 


where 


^’^,o = 





(38) 


After implementing the first unitary operations in the |0)^ part of the 
initial state, we can get normalized state 

Secondly, using the QWD operation once again to transform the |0)i, part 
of initial state into the normalized state We let g = Xlhi 
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and define the second QWD operation as , which can be expressed as a 
L X L matrix. The elements of the matrix satisfy 


T/S _ 






S 12 
ol 


where 




(39) 


(40) 


After implementing the second unitary operations on |0 )l part of initial 
state, we can get normalized state X)Li 

We perform the L level auxiliary qudit |0 )l and K auxiliary qubits |0)^ 
controlled operation Ui on the computer. The processing can be described as 



fc=0 1=1 


(41) 


Then, corresponding to the two times of performing of the QWC opera¬ 
tions, we need perform the QWC operations twice to combine the wave func¬ 
tions. We define the QWC operations IT^and corresponding to the QWD 
operations and respectively. 

We set = (C'®)’!' and perform the the QWC operations 

IT^and on the state and |£), respectively. 

After QWC operations IT^and IT'®, we detect the final wave function when 
the auxiliary system is in state |0)^|0)i,. In the final state, we only focus our 
attention on the terms with the L level auxiliary qudit in state |0 )l and K 
auxiliary qubits in state |0)^. We have the following 


1 

77 




t.—n 


L 


E 




K) 


L 


e.=i 


(42) 

(43) 


It should be noted that the summation parts /k\ and X)7i 

already have been combined with Ui. The initial state is transformed into 

K L 

I'^)|0)'^|0)l ^ /k\Y^aiUm^)^\^)L, (44) 

where Ui corresponds to some (—and i G J. 

It is obvious that s = fg and EjeJ PjVj = Conse¬ 

quently, we have successfully realized the following process: 
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\^)\0)^\0)l^-U\<F)\0)^\0)l. 

S 


(45) 


Finally, the robust form of obvious amplitude amplification procedure of 
[34] enables us to deterministically implement U through amplifying the am¬ 
plitude of —U. The approximation accuracy of U can be quantified by approx- 
s 

imation error e. According to the Chernoff bound, as studied in [^, the query 
complexity is 


log(r/e) \ 
log log(r/e )) ’ 


(46) 


and the error of approximation in each segments satisfy: 


\\U-Ur\\<-. 

r 


(47) 


The total number of gates in the simulation for time t/r in each segment 
is [32 


/ L(n-Hog£)log(T/e) \ 
V loglog(r/e) )■ 


(48) 


In the duality quantum computer description, our method gives a slight im¬ 
provement than [34) . which uses 


/ L(n-b l)log(T/e) \ 
loglog(r/e) J 


(49) 


gates. 

Thus, we have given a standard program in the duality quantum computer 
to realize the simulation methods of Hamiltonians based on linear combina¬ 
tions of unitary operations. The physical picture of our description of the 
algorithm is clear and simple: each of QWD and QWC operations will lead 
to one summation of linear combinations of unitary operations. Our method 
is intuitive and can be easily performed based on the form of time evolution 
operator. 


7 Summary 

In the present paper, we have briefly reviewed the duality quantum computer. 
Quantum wave can be divided and recombined by the QWD and QWC oper¬ 
ations in a duality quantum computer. The divider and combiner operations 
are two crucial elements of operations in duality quantum computing and 
they are realized in a quantum computer by unitary operators. Between the 
dividing and combining operations, different computing gate operations can 
be performed at the different sub-wave paths which is called the duality par¬ 
allelism. It enables us to perform linear combinations of unitary operations in 
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the computation, which is called the duality quantum gates or the generalized 
quantum gates. The duality parallelism may exceed quantum parallelism in 
quantum computer in the precision aspect. 

The duality quantum computer can be perfectly simulated by an ordinary 
quantum computer with n-qubit and an additional qudit, where a qudit labels 
the slits of the duality quantum computer. It has been shown that the dual¬ 
ity quantum computer is able to simulate any linear bounded operator in a 
Hilbert space, and unitary operators are just the extreme points of the set of 
generalized quantum gates. 

Simulating the time evolution of quantum systems or the dynamics of quan¬ 
tum systems is a major potential application of quantum computers. The prop¬ 
erty of duality parallelism enables duality quantum computer to simulate the 
dynamics of quantum systems using linear combinations of unitary operations. 
It is naturally suitable to realize the simulation algorithms of Hamiltonians 
based on multi-product formulas which are usually adopted in classical algo¬ 
rithms. The duality quantum computer can be used as a bridge to transform 
classical algorithms into quantum computing algorithms. We have realized 
both Childs-Wiebe algorithm and Berry-Childs simulation algorithms in a du¬ 
ality quantum computer. We showed that their algorithm can be described 
straightforwardly in a duality quantum computer. Our method is simple and 
has a clearly physical picture. Consequently, it can be more easily realized in 
experiment 

This work was supported by the National Natural Science Foundation of 
China Grant Nos. (III75094 , 9I22I205), the National Basic Research Pro¬ 
gram of China (2011CB92I6002), the Specialized Research Fund for the Doc¬ 
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